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ABSTRACT 

Context. We previously reported the direct detection of a low mass companion at a projected separation of 55±2 AU around the B9 type star k 
Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of 
gas giant planets and brown dwarfs on wide-orbits. 

Aims. We present new angular differential imaging (ADI) images of the system at 2.146 (K s ), 3.776 (L’), 4.052 (NBA. 05) and 4.78 pm (M’) 
obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We 
aim to determine the near-infrared spectral energy distribution (SED) of the companion and use it to characterize the object. 

Methods. We used analysis methods adapted to ADI to extract the companion flux. We compared the photometry of the object to reference 
young/old objects and to a set of seven PHOENIX-based atmospheric models of cool objects accounting for the formation of dust. We used 
evolutionary models to derive mass estimates considering a wide range of plausible initial conditions. Finally, we used dedicated formation 
models to discuss the possible origin of the companion. 

Results. We derive a more accurate J = 15.86 ± 0.21, H = 14.95 + 0.13, K s = 14.32 ± 0.09 mag for k And b. We redetect the companion in all 
our high contrast observations. We confirm previous contrasts obtained at K s and L’ band. We derive NBA. 05 = 13.0 ± 0.2 and M' = 13.3 ± 0.3 
mag and estimate Logi 0 (L/L o ) = -3.76 ± 0.06. Atmospheric models yield T eff = 1900^oo K They do not set constrains on the surface gravity. 
“Hot-start” evolutionary models predict masses of 1 4^ 5 M Jup based on the luminosity and temperature estimates, and considering a conservative 
age range for the system (30^° Myr). “warm-start” evolutionary tracks constrain the mass to M > llM Jup . 

Conclusions. The mass of k Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity 
models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible 
masses of k And b. 
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1. Introduction 

During the last 15 years, radial velocity and transit surveys have 
provided a detailed inventory of the population of giant plan- 
ets within ~3 AU around stars of different masses, ages, and 
metallicities (e.g Lagrange et al. 2009b; Johnson et al. 2011; 
Sousa et al. 2011; Mortier et al. 2012; Bonfils et al. 2013; Sato 
et al. 2013; Nowak et al. 2013, and ref therein). Correlations be- 
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tween planet frequencies and the host- star metallicity (Gonzalez 
1997; Santos et al. 2001; Fischer & Valenti 2005; Sousa et al. 
2011; Mortier et al. 2012), the host-star mass (Lovis & Mayor 
2007; Bowler et al. 2010), and between the heavy-element con- 
tent of gas giants with host star metallicity (Guillot et al. 2006; 
Miller & Fortney 2011) favour the hypothesis of a formation 
by core- accretion (hereafter CA; Pollack et al. 1996; Mordasini 
et al. 2009a, b; Alibert et al. 2011; Mordasini et al. 2012a). Core- 
accretion considers that a core of solids (ice, rock) forms through 
collisions of planetesimals in the protoplanetary disk at a dis- 
tance of a few AU from the central star. Once the core has 
reached a critical mass (Mizuno 1980; Bodenheimer & Pollack 
1986), its gravitational potential causes a rapid capture of the 
surrounding gas which ultimately forms a massive gas envelope. 
Additional migration mechanisms have been proposed (e.g Lin 
& Papaloizou 1986; Alibert et al. 2004, and ref. therein) to ex- 
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Table 1 . Log for the new high-contrast observations of k Andromedae. 


Date 

Instrument 

Band 

DIT 

(s) 

NDIT 

Nexp 

Parallactic angle / start 

O 

Parralactic angle / end 

O 

Remarks 

2012/10/06 

LMIRCam 

M 

0.757 

30 

87 

177.49 

123.54 

Saturated exposures 

2012/10/06 

LMIRCam 

M 

0.029 

4 

36 

165.97 

124.98 

Unsaturated exposures 

2012/10/30 

NIRC2 

L' 

0.3 

70 

30 

160.92 

153.58 


2012/10/30 

NIRC2 

L' 

0.3 

70 

30 

140.35 

134.46 


2012/11/03 

NIRC2 

NB_4.05 

0.3 

100 

10 

163.14 

159.94 


2012/11/03 

NIRC2 

K s 

5 

3 

15 

155.27 

151.92 


2012/11/03 

NIRC2 

NB_4.05 

0.3 

100 

30 

144.95 

137.40 


2012/11/03 

NIRC2 

K s 

5 

3 

15 

133.43 

131.58 



plain the population of giant planets orbiting very close to their 
parent stars. 

Conversely, high-contrast and high-angular resolution imag- 
ing is probing the population of wide-orbit (> 5AU) gaseous 
companions around a variety of young (age <300 Myr) and 
nearby (d < 150 pc) stars, ranging from M dwarfs to early- 
type (F to B) stars (Neuhauser et al. 2002; Chauvin et al. 2003; 
Lowrance et al. 2005; Masciadri et al. 2005; Biller et al. 2007; 
Kasper et al. 2007; Lafreniere et al. 2007a; Chauvin et al. 2010; 
Janson et al. 2011; Delorme et al. 2012; Bowler et al. 2012; 
Vigan et al. 2012; Rameau et al. 2013b; Nielsen et al. 2013). 
The majority of planetary mass companions have been discov- 
ered along an extended range of wide orbits (15 AU to several 
hundreds of AU) around low mass (MGK) stars (e.g Chauvin 
et al. 2004; Todorov et al. 2010; Lafreniere et al. 2010). The 
high-mass ratio with their host and the high separations makes 
the fragmentation of pre-stellar cores during collapse (e.g. Bate 
2012) a candidate for the formation of these wide systems. 

Low mass (< 15 M Jup ) gaseous companions discovered 
more recently at moderate separations (< 100 AU) around 
the massive stars HR 8799 (Marois et al. 2008, 2010) and p 
Pictoris (Lagrange et al. 2009a, 2010) might represent a pre- 
viously unexplored population of gaseous companions (Vigan 
et al. 2012; Rameau et al. 2013b; Nielsen et al. 2013). The 
extended debris disks identified around these stars, shaped by 
the companions (Su et al. 2009; Lagrange et al. 2012a), suggest 
these systems emerged from a primordial gaseous disk. This 
picture is reinforced by recent resolved images of transition 
disks around Herbig stars (Andrews et al. 2011; Rameau et al. 
2012) with cavities which extend beyond the separations of the 
aforementioned companions, and might have been carved by 
planets (Casassus et al. 2013). Nevertheless, the subsequent 
formation pathway remains unclear. Core accretion (CA) could 
eventually explain the properties of directly imaged planets with 
the narrowest orbits (HR 8799 e and d, ft Pictoris b; Kennedy 
& Kenyon 2008; Mordasini et al. 2009a; Rafikov 2011). But 
associated CA formation timescales become too long compared 
to the mean lifetime of primordial disks and require higher 
disk surface density for an in- situ formation at more than ~15 
AU (Boley 2009; Dodson-Robinson et al. 2009). A revision 
of the way solids are accreted (Ormel & Klahr 2010), or 
additional outward migration mechanisms must be considered 
(Crida et al. 2009a; Kley & Nelson 2012, and ref therein) to 
explain planets found at larger radii if formed initially by CA. 
Gravitational instability within disks (G.I.; Cameron 1978) 
has been considered as an alternative mechanism for these 
objects (Boss 2011) and can also be associated with migration 
(e.g. Zhu et al. 2012) and ejection (e.g. Vorobyov 2013). Here, 
protostellar disks develop global instabilities (if cool enough) 
and fragment into bound clumps that contract to form giant 
planets. This mechanism operates on much shorter timescales 


than CA (a few orbital periods). However, recent surveys 
suggest it might not dominate at wide (> 30 AU) separations 
(Janson et al. 2011, 2012; Rameau et al. 2013b). To conclude, 
it is not clear how these 8-78 AU companions relate to the 
more distant low-mass brown-dwarfs companions found around 
other massive (1.35-2.5) young (age < 150 Myr) stars (HR 
7329 B, HIP 78530 B, HD 1160 B, HR 6037 BaBb, HD 23514 
B; Lowrance et al. 2000; Lafreniere et al. 2011; Nielsen et al. 
2012; Huelamo et al. 2010; Nielsen et al. 2013; Rodriguez et al. 
2012 ). 

Direct imaging offers the possibility to collect multiple-band 
photometry and spectra emitted by the photospheric layers of the 
companions in the near-infrared (1-5 /mi; McElwain et al. 2007; 
Currie et al. 2011; Bonnefoy et al. 2010, 2013b). These data can 
provide a stringent characterization of the chemical (composi- 
tion) and physical properties (mass, radius, effective tempera- 
ture) of the sources, which are at the basis of our understanding 
of their formation processes (Bonnefoy et al. 2013a; Konopacky 
et al. 2013). They can also give glimpses of the physics and 
chemistry at play in the cool and complex atmospheres of the 
sources (e.g. Currie et al. 2011; Skemer et al. 2012; Bonnefoy 
et al. 2013b). The peculiar near-infrared spectro-photometric 
properties of the companions are now better understood as a con- 
sequence of the low temperature and surface gravity atmosphere, 
which in some cases can lead to the formation of thick layers 
of dust and/or trigger non-equilibrium chemistry (Janson et al. 
2010; Skemer et al. 2011; Barman et al. 201 la, b; Madhusudhan 
et al. 2011; Currie et al. 2011; Skemer et al. 2012; Marley et al. 
2012; Faherty et al. 2013). 

In the course of the Strategic Explorations of Exoplanets 
and Disks with Subaru (SEEDS, Tamura 2009), we identified 
in early-2012 a low-mass companion at a projected separation 
of 1” around the massive (2.5 ±0.1 M 0 ) B9 IVn (Garrison & 
Gray 1994) and nearby (51.6 ± 0.5 pc, van Leeuwen 2007) star 
k Andromedae (hereafter k And, Carson et al. 2013). The host 
star kinematics make it a high probability member (Zuckerman 
et al. 2011; Malo et al. 2013; Carson et al. 2013) of the ~30 Myr 
old Columba moving group (Torres et al. 2008). At these ages, 
we estimated k And b has a mass of 12.8^^ M Jup , making it the 
first planetary/brown-dwarf companion directly imaged around 
such a massive star. 

However, the mass estimate of k And b relied on the near- 
infrared photometry of the star, derived from 2MASS images 
(Cutri et al. 2003) where k And A is saturated. More importantly, 
it was derived from predictions of the so-called “hot-start” evo- 
lutionary models, which assume that the object starts its evolu- 
tion following a spherical collapse from an arbitrary large ini- 
tial radius (Burrows et al. 1997; Chabrier et al. 2000; Baraffe 
et al. 2003). The alternative “cold- start” models (Marley et al. 
2007; Spiegel & Burrows 2012; Mordasini et al. 2012b) hypo- 
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thetize that the gaseous material accreted onto planet embryo 
passes through a super-critical accretion shock, and looses all its 
gravitational energy, therefore leading to objects with low ini- 
tial entropies As a consequence, these models predict 

lower temperatures and lower luminosities at early ages than 
“hot- start” evolutionary models for a given objet. The more re- 
cent “warm-start” models (Spiegel & Burrows 2012; Marleau 
& Cumming 2013) generalise the previous cases by exploring 
the impact of initial conditions on the cooling curves through 
the choice of the initial entropy (Si n i t ) of the object. With these 
models, joint constraints on the mass and initial entropy of com- 
panions can be derived from a brightness and temperature mea- 
surement. 

We present new high-contrast near-infrared images of k And 
b obtained using NIRC2 at the W.M. Keck Observatory and 
LMIRCam on the Large Binocular Telescope Interferometer 
(LBTI) at 2.146 (K s ), 3.776 (L’), 4.052 (NB_4.05) and 4.78 pm 
(M’). They confirm and complement the current set of photomet- 
ric data of the companion. We also present additional unresolved 
observations of the system in the near-infrared. Unresolved ob- 
servations provide accurate photometry of the primary star and, 
as a consequence, of the companion from 1 to 2.5 pm, at wave- 
lengths where the effect of atmospheric dust can be studied. 
Altogether, they enable us to refine the companion properties 
and discuss its formation mechanism. 

This paper is organized as follows: we describe in Section 2 
the observations and the related analysis of the data; we present 
in section 3 our main results. Section 3 is split into four subsec- 
tions. We first rederive a more conservative age estimate for the 
system in Subsection 3.1. We compare the photometry of k And 
b to empirical reference objects in Subsection 3.2, and to atmo- 
spheric models in Subsection 3.3. We give new mass estimates 
based on “hot”, “cold”, and “warm-start” evolutionary models in 
Subsection 3.4. We discuss the properties of k And b and review 
the possible formation scenarios in section 4. 

2. Observation and data reduction 

2.1. Seeing-limited observations 

The accuracy of k And b near infrared photometry (J=16.3±0.3, 
H=15.2 ± 0.2, K s = 14.6 ± 0.4) reported in Carson et al. (2013) 
was predominantly limited by the accuracy of the 2MASS pho- 
tometry of the star (the only near-infrared photometry publicly 
available at that time; Table 2). The low accuracy arises from the 
strong saturation of the star in the 2MASS images 1 . 

We obtained infrared photometry of k And A in the J, H, and 
K s filters using the MIMIR instrument (Clemens et al. 2007), 
mounted on the 1.8m Perkins telescope at Lowell Observatory, 
on 2012 October 30. At the f/5 focus, the 1024 x 1024 Aladdin 
array covered a field of view (FOV) of 10 x 10 arcmin with a 
pixel scale of Off 519 pixel -1 . Standard bias, dark, flat field, and 
sky corrections were applied to all data using calibration data 
obtained during the night. Five 0.1s exposures were obtained at 
each of six dither positions separated by -60 arcsec in every fil- 
ter. Due to the brightness of our science target, we defocus sed 
the telescope until stars took the appearance of -60 arcsec wide 
donuts. Our observations of k And A were bracketed by obser- 
vations of HR 8962, using the same exposure time and focus 
settings, to enable us to determine the photometric zeropoint. 

Photometry was determined using apertures with a radius of 
70 pixels and a 10 pixel wide background sky annulus starting 

1 The saturation can be seen directly on the 2MASS images: 
http://irsa.ipac.caltech.edu/cgi-bin/2MASS/IM/nph-im_pos 


Table 2. Near-infrared photometry of k And A 


Epoch 

J 

H 

K s 

Ref. 

29/10/1998 

4.62 ± 0.27 

4.60 ± 0.22 

4.57 ± 0.36 

1 

30/10/2012 

4.26 ± 0.04 

4.31+0.05 

4.32 ± 0.05 

2 


Notes. References: [1] - 2MASS / Cutri et al. (2003), [2] - this work 


at a radius of 80 pixels from each centroid position. Our final 
photometry was computed from 10 data frames (2 dither posi- 
tions) in which our large apertures did not intersect with known 
regions of bad pixel clusters in the array. We determined the ze- 
ropoint for our observations by comparing our measured pho- 
tometry for HR 8962 against published 2MASS values (Cutri 
et al. 2003). HR 8962 is known to be a visual binary (see e.g. 
Morlet et al. 2002; Hartkopf & Mason 2009) whose separa- 
tion of 07 5 is within both our aperture definition and the aper- 
ture definition of 2MASS, and thus should not affect our zero- 
point determination. The photometry we determined for 2MASS 
J23373296+4423122, another bright object in the FOV of our 
HR 8962 data, was within 2 -cr of published 2MASS photometry 
for the source, thereby confirming our zeropoint determination. 

Aperture photometry on the source and of the photo- 
metric reference yields J = 4.26 ± 0.04, H = 4.31 + 0.05, and 
K s = 4.32 ± 0.05 mag for k And A (see also Table 2). We find 
negligible photometric shifts (< 0.004 mag) between MIMIR 
and HiCIAO photometric systems using the corresponding filter 
pass-bands, a flux-calibrated spectrum of Vega (Bohlin 2007), 
and of a B9 star from the Pickles (1998) library. This enables us 
to revise the original photometry of k And b (Carson et al. 2013) 
to J = 15.86 ± 0.21, H = 14.95 ± 0.13, K s = 14.32 ± 0.09 mag. 

2.2. High-resolution spectroscopy 

We obtained a R -31500 optical (-3600-10000A) spectrum 
of k And A on UT 2012 October 24 with the ARC Echelle 
Spectrograph (ARCES, Wang et al. 2003) mounted on the 
Apache Point Observatory (APO) 3.5 m telescope. The spec- 
trum was obtained using the default 176 x 372 slit and an ex- 
posure time of 45 seconds. A ThAr lamp exposure was obtained 
after the integration to facilitate accurate wavelength calibration. 
The data were reduced using standard IRAF techniques. We used 
this spectrum to derive new estimates of the surface gravity and 
effective temperature of the star. Our results and method are de- 
tailed in Section 3 and Appendix B. 

2.3. High-contrast observations 
2.3.1. Keck/NIRC2 

We observed k Andromedae on October 30, 2012 and November 
3, 2012 with the NIRC2 camera fed by the adaptive optics 
system of Keck II (van Dam et al. 2004). We used the K s 
(A c = 2.16//m) and L' (A c = 3.78/mi) broadband filters and the 
Brer (A c = 4.05/rm) narrow band filter (hereafter NB_4.05, Table 
1). All data were taken in the narrow camera mode (9.952 
mas pixel -1 , Yelda et al. 2010) in either multi-correlated dou- 
ble sampling (K s ) or correlated double sampling (L', NB_4.05). 
Observing conditions both nights were photometric with above 
average seeing (FWHM natura i ~ 074—075) and average AO per- 
formance. In all cases, we used the “large hex” pupil plane mask. 

The K s data were taken through the partially transmissive 
07 6 diameter coronagraphic mask in coadded 15 second expo- 
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Fig. 1 . Reduced Keck/NIRC2 images at K s /2.16 pm (left), L'/3.78 pm (middle), and NB_4.05 (right) of k And showing the detection 
of the companion k And b. We adjust the image color scale such that the white region identifying k And b (pointed to by an arrow) 
roughly covers an area equal to the image FWHM. The cross identifies the stellar centroid position in the K s data and is at the same 
point in all three images. 


sures. The L' data were taken in coadded 20 second exposures 
with the primary star’s point- spread function (PSF) core satu- 
rated. We dithered the star off of the detector to obtain sky frames 
in the middle of our observing sequence. For the NB_4.05 data, 
the PSF core was unsaturated, and we took three sequences of 
10 science frames followed by 5 sky frames (t exp = 30 s). 

All data were taken in “vertical angle” or angular differential 
imaging mode (Marois et al. 2006). On both nights, we observed 
k And immediately after unsaturated observations of HR 8799, 
which is at a similar right ascension, beginning at an hour angle 
of ~ 0.5. These unsaturated observations were later used to de- 
rive the photometry of k And b (see below). Still, over the course 
of our sequence, we achieved ~ 24-26° of field rotation, or « 5- 
10 T/D at the angular separation of k And b. 

Basic image processing followed standard steps previously 
used to process NIRC2 data (Currie et al. 2012b,c). Briefly, 
for the K s data, we employ standard dark subtraction and flat- 
fielding corrections, identify and interpolate over hot/cold pixels. 
For the thermal IR data, we subtracted off a median-combined 
sky image comprised of sky frames taken closest in time to the 
science frames of interest. For all data sets, we applied the distor- 
tion correction from Yelda et al. (2010). After copying each im- 
age into a larger blank image, we performed image registration 
by finding directly the stellar centroid position (K s and NB_4.05) 
or estimating it by cross-correlating the first image in our se- 
quence by a 180° rotation of itself (L'). We then determined the 
relative offsets for other images in the sequence by solving for 
the peak in the cross-correlation function between the first (ref- 
erence) image and all subsequent images. 

To extract a detection of k And b, we used the A-LOCI 
pipeline described in Currie et al. (2012a). Because k And b 
is at a wide separation (r 1705) and previous data obtained 
at similar wavelengths yielded very high SNR detections, we 
adopt conservative algorithm settings that minimally bias the 
planet flux. This included a large rotation gap (£ = 0.8-1. 5), a 
high cross-correlation cutoff for the contrast-limited K s and L' 
data ( corrum =0.9-0.95) and, a lower cutoff for the background- 
limited NB_4.05 data (corru m = 0.2), and a moving pixel mask 
(Currie et al. 2012a), yielding planet throughputs (estimated by 
implanting synthetic point sources into registered images) rang- 
ing between 0.92 and 1. 

Figure 1 shows the final Keck/NIRC2 images, where we de- 
tect k And b at a SNR of 30, 20, and 7 at K s , L', and NB_4.05. 
Exterior to our inner working angle of O'/ 3, we do not detect any 
other point sources, as was the case with previous H, K s , and L' 
data (Carson et al. 2013). To measure the brightness of k And b, 


we perform aperture photometry with a diameter roughly equal 
to the image FWHM (-5-10 pixels), correcting for the very mi- 
nor throughput loss induced by our processing. To flux calibrate 
k And b in niK s , we use the derived coronagraph spot extinction 
from Currie et al. (2012b) of 6.91 ± 0.15 mag. This yielded a 
parent/companion contrast of Am = 10.04 ± 0.15 mags. Using 
our MIMIR photometry, this translates into a companion bright- 
ness of m Ks = 14.36 ± 0.15. For the L', we used the unsaturated 
images of HR 8799 A (m NB _ 4 .o 5 - m L / = 5.220 ± 0.018; Marois 
et al. 2008) to flux-calibrate k And b, deriving m L / = 13.13 ± 
0.07. Finally, for the NB_4.05 data, both HR 8799 A and k And 
A were unsaturated. Assuming m L > ^ m NB _ 4 .o 5 , we then derive 
mNB_ 4.05 = m L / = 4.32 ± 0.05 for k And A and m NB _ 4 . 0 5 = 13.0 
± 0.2 for k And b. 

We measure a separation of p = 1.029 ± 0.005” and a po- 
sition angle of 6 = 55.3 ± 0.3° in the K s band images. This 
value is in agreement with the astrometry derived from IRCS 
(p = 1.044”, 6 = 55.2) and HiCIAO data presented by Carson 
et al. (2013). We nevertheless refrained from making an updated 
proper motion analysis, due to possible systematic offsets on p 
and 6 introduced by the instrument change. 

2.3.2. LBTI/LMIRCam 

We observed k And with the LMIRCam near-infrared camera 
(Hinz et al. 2008; Skrutskie et al. 2010) at the LBTI on October 
10, 2012. The LBTI was operated in single- aperture mode in or- 
der to avoid extra-overheads associated with the alignment of the 
telescope beams. The telescope+instrument do not have a dero- 
tator. Therefore, it automatically operates in a mode that enables 
passive angular differential imaging (Marois et al. 2006). We ob- 
tained 151 frames consisting of 30x0.758s coadded exposures 
each with a M’-band filter ( A c = 4.78pm, FWHM = 0.37pm). 
The integration time was chosen in order to saturate the core of 
the stellar point spread function (PSF) to a radius of 110 mas. 
The lh45 min spent on the target produced a field rotation of 
53.9°. The source was nodded in the instrument field of view 
(1.8” nod) every 3 mins in order to properly remove the back- 
ground contribution. We recorded four frames corresponding to 
eight 0.029s coadded exposures following each telescope nod. 
These unsaturated exposures monitored the evolution of the PSF 
during the observing sequence, and were used to derive the con- 
trast ratio of the system components. 

Data were reduced with the MPIA-LBTI angular differen- 
tial imaging pipeline. The pipeline carried out all the basic cos- 
metics steps (removal of detector stripes, sky subtraction, bad 
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Fig. 2. LBTI/LMIRCam high-contrast image of the k And sys- 
tem at M band. The companion is recovered. 


pixels interpolation) following the selection of pairs of nodded 
exposures. The position of the source in the resulting images 
was found using a bidimensional Moffat fitting function. Images 
were recentered and placed in a master cube of frames using 
subpixel shifts. We selected 87 frames with high-Strehl ratio or 
high-throughput to build a final master cube. We applied six 
algorithms to estimate and remove the flux distribution of the 
star in each input frame of the master cube. We removed a ra- 
dial profile in each image, derotated them in order to re-align 
them with the North using the parallactic angle at the time of 
the exposure (non- ADI, or NADI), and median-combined these 
frames to produce a final residual image. Alternatively, we built 
and removed the PSF of the star taking the median of all in- 
put frames contained in the cube (classical ADI, or CADI). We 
also applied the RADI (Marois et al. 2006; Chauvin et al. 2012), 
LOCI (Lafreniere et al. 2007b), and a custom algorithm based 
on Principal-Component Analysis (Zimmerman et al., in prep; 
see also Soummer et al. 2012; Amara & Quanz 2012). 

We also developped and applied an extension of the RADI 
algorithm, called “correlated radial ADI” (or CODI), in order 
to optimise the reduction of these data. Each cube image i was 
first divided into concentric annuli. For each annuli ANi , we se- 
lected N reference annuli AN re f k= \ _ N taken from other images 
contained in the cube for which the companion has moved by 
nx FWHM. We computed and sorted the correlation factors be- 
tween the N reference annuli and the processed annuli. We se- 
lected the a (e N) most correlated reference annuli AN re f k=l _ a 
and retrieved the coefficients of the linear relation between pix- 
els of the reference annuli and of the processed annuli ANi = 
ak,i x AN re f k + bk,i (with a kj and bk,i e R. derived from a linear 
regression and a matrix inversion). This way, we were able to ac- 
count simultaneously for any variation of the atmospheric trans- 
mission and of a piston flux offset (e.g. residuals from the back- 
ground subtraction, etc) between the reference and processed an- 
nuli. The a model annuli were median-combined to create a final 
model AN mo dei of AN t . The operation was repeated for each con- 


centric annuli and each cube image. We considered a separation 
criterium of 1.5 FWHM at the separation of the companion for 
the CADI, RADI, CODI, LOCI, and PCA-based analyses. 

k And b is detected in each individual input frame. It is then 
naturally retrieved in all post-processed frames with a signal to 
noise from 20 to 38 reached with the RADI and CADI algo- 
rithm, respectively (see Figure 2). We integrated the flux of k 
And b over an aperture of 16 pixels in radius (1.5 FWHM) and 
used the non-saturated exposures to derive the contrast ratio be- 
tween the star and its companion. Values were corrected from 
the inevitable flux losses associated to angular differential imag- 
ing using 3 fake planets injected at position angles of -67°, 115°, 
and 233°. 

All algorithms converge to a contrast of AM' = 8.9 ± 0.3 
mag for k And b. The variation of the PSF in unsaturated 
exposures dominates the final error budget. We estimate that 
k And A has a M-band magnitude of 4.4 mag consider- 
ing the mean K s - M colors of B9 (IV- V) stars of the van 
der Bliek et al. (1996) and IRTF (Leggett et al. 2003) cata- 
logues (K s - M = -0.023 ± 0.030). We also find a similar color 
(K s - M = -0.023 ± 0.006) using ATLAS9 models (Castelli 
et al. 1997) taken at the temperature of the star (see section 3.1 
and Appendix B), the K s , and M’ band filter transmission curves. 
These colors can be used as a reasonable guess for k And A be- 
cause of the lack of excess emission for this star at these wave- 
lengths (see Section 3.1). Therefore, we derive M’ = 13.3 ± 0.3 
mag k And b. 

We did not get observations of an astrometric field necessary 
to derive a reliable astrometry from these data. 

3. Analysis and Results 

3.1. Re-evaluation of the system age 

The mass of k And A and its companion, and insights into 
the companion’s origin, are tied to the correct determination 
of the system age. For this reason, we re-investigated the 
various age indicators available for k And A based on our new 
measurements of the star and material found in the literature. 

Carson et al. (2013) suggest a possible age range of 20 - 120 
Myr, and a nominal age of ~30 Myr, for the system based on 
the kinematics and color-magnitude diagram (CMD) position of 
k And A. Their kinematic study was based on the Malo et al. 
(2013) online tool 2 , which computes probabilities of member- 
ship in the TW Hydrae, Tucana-Horologium, Columba, Carina, 
and Argus associations, the ft Pictoris and AB Doradus mov- 
ing groups, and the field population using Bayesian methods. 
We found a 95.6%, 0.7%, and 3.8% chance that k And belongs 
to Columba, ft Pictoris, and the field, respectively. The star has 
a 0% probability to belong to the remaining groups. The same 
tool applied to HR8799 yields 98.1%, 0.7%, and 1.1% prob- 
ability to Columba, p Pictoris, and the field. Recently, Baines 
et al. (2012) put independent constraints on the age of HR8799 
which are consistent with the high probability of Columba mem- 
bership derived using this tool. One might argue that the Malo 
et al. (2013) Bayesian analysis tool assumes k And is a mem- 
ber of Columba in the priors of the calculation. This will artifi- 
cially inflate the probability of group membership since the star’s 
kinematics are partly being used to define those of the group. 
An independent analysis of the k And’s UVW space velocities 
shows that they are consistent with other proposed group mem- 
bers with measured parallaxes at the <2cr level, k And’s X and 

2 http://www.astro.umontreal.ca/~malo/banyan.php 
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Fig. 3. Position of k And A in a color-magnitude diagram (CMD) 
compared to other A and B-type member of the Pleiades, the IC 
2391 cluster, Scorpius-Centaurus (ScoCen), UrsaMajoris group, 
and other young moving groups in the solar neighborhood (see 
legend). The star’s position in the diagram is consistent with 
some members of all the plotted samples except the Ursa Majoris 
moving group (UMa), which is clearly older and more evolved. 
The placement in the CMD suggests that k And is not much older 
than the upper limit of the Pleiades, 150 Myr. The dashed-dotted 
lines are Ekstrom et al. (2012) tracks that include rotation. These 
models predict an age of <250 Myr for k And, but they overes- 
timate the ages of B-type cluster members. The dotted lines are 
5th order polynomial fits to the Pleiades and IC 2391 sequences 
and are shown to highlight them for clarity. 


Z Galactic distances are also completely consistent with those 
of the previously mentioned members with well constrained dis- 
tances. The star’s Galactic Y distance (46.5 pc) however, falls 
above the mean value for Columba’s bona fide members (-26.3 
pc; using Table 3 of Malo et al. 2013, and removing k And from 
the sample). While discrepant from the group members origi- 
nally proposed in Torres et al. (2008), this is consistent with 
other proposed Columba members from Zuckerman et al. (201 1) 
that lie at northern declinations. Additionally, the dispersion in Y 
values for proposed members is large (26.5 pc) and might be un- 
derestimated since few surveys have searched for new members 
in the north. Thus, there are several lines of evidence that support 
k And’s kinematic membership in the proposed Columba asso- 
ciation. A full kinematic traceback study of proposed Columba 
group members may shed more light on k And’s reliability as a 
member and the past history of the association as a whole (e.g. 
Ortega et al. 2007; Makarov 2007; Weinberger et al. 2013). Since 
consistent kinematics are a necessary but not sufficient criteria 
for moving group membership, we investigate the age of k And 
using several different methods that are independent of its kine- 
matics. 

The existence of excess IR flux in the spectral energy distri- 
bution (SED) of k And A is indication of remnant, circumstel- 
lar material (a debris disc) and may place constraints on its age 


Fig. 4. Position of k And A in a temperature-gravity diagram 
compared to the models of Ekstrom et al. (2012) for rotating 
(solid lines; V/V cnr =0.4 when stars reach the zero-age main- 
sequence) and non-rotating (dashed limes). We also show the 
same parameters for similarly typed members of young clusters 
(see legend). The models do not reproduce the ages of the clus- 
ter members and there is significant scatter in the polar gravity 
estimates within a given cluster. We interpret this as indication 
that these comparisons do no place a meaningful constraint on 
the age of k And. 


(see discussions in Wyatt 2008). To construct k And A’s SED, we 
used the Virtual Observatory (VO) SED Analyzer (VOSA, Bayo 
et al. 2008) 3 . The tool allows the user to use both publicly avail- 
able and user provided data to construct the SED of a source 
and fit it with their choice of model. We queried photometric 
catalogs available through the VO to compile a complete SED 
of k And A. We recovered data from -0.13 to 100 //m from the 
following sources: the International Ultraviolet Explorer (IUE, 
Boggess et al. 1978), A catalogue of compiled UBV photom- 
etry (Mermilliod & Mermilliod 1994), the Tycho-2 catalogue 
(H0g et al. 2000), the 2MASS All-Sky Point Source Catalog 
(Cutri et al. 2003), the WISE All-Sky Data Release (Cutri & 
et al. 2012; Wright et al. 2010), the AKARI/IRC mid-IR all-sky 
Survey (Ishihara et al. 2010), and the IRAS Catalog of Point 
Sources, Version 2.0 (Helou & Walker 1988). Not all of the 
recovered data were useful for constructing the SED. As men- 
tioned in Section 2.1, the 2MASS photometry of k And A is sat- 
urated, we therefore replaced it with our new determinations of 
the JHK S photometry. The WISE W 1 and W2 data were also be- 
yond the saturation limit listed in the Explanatory Supplement 4 , 
so we did not include these points in the construction of the SED. 
The IRAS 25, 60, and 100 pm data are listed as upper limits, 
thus, they were also not used. To supplement the SED, we also 
checked for additional mid-IR data in the S pitzer and Herschel 
archives (not queried by the VO). The star was not observed by 
either telescope. We therefore used the VOSA to perform a x 1 fit 


3 http://svo2.cab.inta-csic.es/theory/vosa/ 

4 http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/ 
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of an ATLAS9 model (Castelli et al. 1997) to the remaining reli- 
able photometric data. The SED fit reveals no significant excess 
above the expected photospheric flux in an 11000 K model out 
to 22 /jm. Unfortunately, this imposes no constraint on the age 
of k And A. Wyatt (2008) shows the evolution of 24 and 70 fim 
excess in A- type stars (a reasonable proxy for k And’s B9 type) 
with ages up to 800 Myr (his Figure 6). The fraction of stars with 
measured excess is a function of age, however more than half of 
the stars observed at 24 /urn have no detectable excess. Even at 
70 fim, not all of the young targets presented in Wyatt (2008) ex- 
hibit excesses. Thus, longer wavelength observations (>22 fim) 
may yet reveal IR excess in the SED of k And A, but the absence 
of excess in the current data does not provide useful information 
regarding its age. 

k And is the earliest- type (B9) proposed member of the 
Columba association. The ages of later-type stars in the asso- 
ciation are well constrained by a combination of CMD positions 
and lithium depletion. These diagnostics indicate an association 
age of -30 Myr Torres et al. (2008). Age determinations like 
lithium depletion are not applicable to young, early-type stars. 
Therefore, to constrain k And’s age independently of its kinemat- 
ics, we rely on HR diagrams and model comparisons. We first 
place k And A in a M v vs. B-V CMD to critically compare its 
position to those of early-type members of several young, open 
clusters with well defined ages (see Figure 3 and Appendix A). 
k And’s CMD position is consistent with similarly typed mem- 
bers of the Scorpius-Centaurus subgroups (ScoCen, 11-17 Myr, 
Chen et al. 2012) and the Pleiades (130 ± 20 Myr, Barrado y 
Navascues et al. 2004). There are no members of the IC 2391 
cluster (30-50 Myr, Barrado y Navascues et al. 2004; De Silva 
et al. 2013) with CMD placement in the immediate vicinity of 
k And, however the approximate shape of the cluster sequence 
would place them close, k And A is clearly younger than the 
-500 Myr Ursa Majoris moving group (King et al. 2003), where 
the earliest type stars are much more evolved. The other pro- 
posed members of young moving groups (Malo et al. 2013) oc- 
cupy small regions of color space and exhibit significant scatter 
in the CMD, therefore comparison to them cannot effectively 
constrain the star’s age. The sample of ScoCen stars is also scat- 
tered. This may be partly attributed to large distance uncertain- 
ties. This aspect of the comparison sample may skew the in- 
terpretation of k And’s age to younger values. Thus, while the 
CMD position of k And is consistent with populations of stars 
at several different ages, we infer that the star’s placement is 
most consistent with the Pleiades. We also compare the CMD 
positions’ of the stars to the rotating models (V/V cn Y=0.4 a t the 
zero-age main- sequence) of Ekstrom et al. (2012). The models 
reproduce the distribution of A-type stars on the main- sequence, 
but, they deviate from the known ages of stars that are beginning 
to evolve to redder colors. For example, B-type IC 2391 and 
Pleiades members are approximately coincident with the 100- 
160 Myr and 160-250 Myr isochrones, respectively. The age of 
k And is predicted to be <250 Myr. However, since member ages 
of well defined clusters are overestimated by the models, we do 
not use them to place limits on the age of k And. Rather, we use 
the comparison to the empirical Pleiades sequence to suggest an 
upper age limit of -150 Myr for the system. 

We also made new estimates of the atmospheric parameters 
of k And A (see the details in Appendix B) and compared these 
values to the tracks of Ekstrom et al. (2012) in both the rotating 
and non-rotating cases. Figure 4 shows our measured polar sur- 
face gravity and effective temperature for the star (see Table 6). 
The polar surface gravity is a correction to the measured surface 
gravity that accounts for the rapid rotation of k And (vsin/=130- 


190 km s“\ Huang & Gies 2006; Glebocki et al. 2000; Abt et al. 
2002; Fitzpatrick & Massa 2005). We also plot the same param- 
eters for samples of mid to late-B type stars in other young clus- 
ters with ages -10-60 Myr (Huang et al. 2010; Marsh Boyer 
et al. 2012). The models predict an age >250 Myr for k And. 
However, they systematically overestimate the ages of all but a 
few of the cluster members. There is also large scatter in the esti- 
mated polar gravity of stars within the same cluster. We interpret 
these features as indication that placement in this kind of dia- 
gram does not provide a meaningful constraint on the age of k 
And. The figure also illustrates the inherent difficulty in accu- 
rately measuring atmospheric parameters of young, early type 
stars. 

The age suggested by the empirical CMD is moderately 
older than that suggested by Columba moving group member- 
ship (30 ^q Myr, Marois et al. 2010) and used by Carson et al. 
(2013). A possible explanation that could reconcile the age es- 
timated following these two approaches is unresolved binarity 
of the star. Future monitoring of the radial velocity or observa- 
tions at higher angular resolution (i.e. sparse aperture masking) 
may be able to constrain this hypothesis. Our high resolution 
spectrum does not resolve k And as a tight binary. However, if 
it were a binary, and the period were long enough, the radial 
velocity amplitude between resolved component lines could be 
overshadowed by rotational broadening in the spectrum. Another 
possible explanation is large intrinsic scatter in the observed pa- 
rameters of early-type stars due to the effects of inclination and 
rapid rotation. When a star rotates rapidly, centrifugal force leads 
to a deformation of the photosphere, and thus a surface grav- 
ity gradient, between the equator and poles. This gives rise to a 
temperature gradient across the observed stellar surface where 
the equator is cooler than the pole. This is known as “gravity 
darkening” (von Zeipel 1924; Huang & Gies 2006, and refer- 
ences therein). Consequently, the measured temperature, surface 
gravity, and other observables are dependent on the viewing an- 
gle of the rotation axis with respect to the observer. The end re- 
sult of these intrinsic effects in early-type stars is that positions 
in diagrams comparing observed parameters (e.g. temperature- 
gravity, temperature-luminosity, color-magnitude) are degener- 
ate in mass, age, and rotation. 

Evidence of these effects has been directly observed in inter- 
ferometric and spectroscopic measurements of the B7V/B8IV 
star Regulus ( a Feo, McAlister et al. 2005; Che et al. 2011). 
The authors find the star is observed edge on and has vsin / ^320 
km s -1 . This combination of rotation and inclination angle lead 
to an observed >3000 K temperature gradient between the pole 
and equator. As a result, the true luminosity of the star is larger 
than that estimated from photometry. Once the effects of rotation 
and inclination are taken into account, the HR diagram age of 
Regulus is reduced by nearly 100 Myr (Che et al. 2011). These 
results may help to reconcile the previously estimated age differ- 
ence between Regulus and its -1767 K2Y companion. Gerbaldi 
et al. (2001) estimate ages for the primary and secondary of 
150 Myr and 30-50 Myr respectively from model luminosity- 
temperature diagrams. Gies et al. (2008) also discovered a low- 
mass, short period (-40 days) companion to Regulus using long 
term radial velocity monitoring. They propose the companion is 
either a white-dwarf or an M dwarf. If the companion is a white 
dwarf, correct age determination of the system becomes more 
difficult because the evolution of the progenitor must be consid- 
ered (Rappaport et al. 2009). These deep investigations of the 
Regulus system may shed light on the observed scatter in sam- 
ples of known age in Figure 3 and Figure 4. The estimated vsin/ 
of k And is not as large as that of Regulus, but the previous ex- 


7 


M. Bonnefoy et al.: Characterization of the gaseous companion k Andromedae b 


: ▼ M0-M4 
3-0 ^ v M5-M9 

> L0-L4 

> L5-L9 


2.5; | 

; : 

2.0 1 
Li 

* 1.5; 
1 . 0 ; 
0.5; 

0.0 L 


I T0-T5 
] T5-Y0 



HD1160B A HR8799b* 
HD1 160 C □ HR8799 c ♦ 

HR8799 d A 
HR8799 e a 
1RXS1609 b • 
AB Pic b a 
HIP78530 B ♦ 



2M1207 A ► 
2M1207 b i 
GSC06214 b > 
2M2249 A ◄ 
2M2249 B x 
(3 Pic b - 
2M0103(AB)b o 

ffl 



- 0.1 


-0.5 0.0 

L’ - M’ 


0.5 


1 

J-K s 


-1.0 -0.5 0.0 


0.5 
H - K 


1.0 1.5 


2.0 


Fig. 5. k And b (magenta dot) in color-color diagrams compared to known MLTY field dwarfs and brown-dwarfs, and young directly 
imaged companions. The photometry of young companions is referenced in Bonnefoy et al. (2013a,b) (and ref. therein), Nielsen 
et al. (2012), and Delorme et al. (2013). The photometry of MLTY field dwarfs is taken from Leggett et al. (2010, 2013). 


ample highlights that rotation/inclination induced effects on the 
measured physical parameters may be significant and change the 
age interpretation. Further examples of interferometric studies of 
rotating, early-type stars and discussions of how new data led to 
revised understandings of long studied stars can be found in van 
Belle (2012, and references therein). 

Future observations of k And and its companion will pro- 
vide the means to refine the system age. The proximity and lu- 
minosity of k And make it amenable to a full interferometric 
analysis where the measured radius, oblateness, and inclination 
can place strong constraints on its evolutionary status (van Belle 

2012) . Additionally, direct spectroscopic observations of k And 
b should allow for independent age constraints via gravity sen- 
sitive features. In the absence of further observations, we use in 
the subsequent analyses: 1/ An age for the k And system based 
on the proposed kinematic membership to the Columba associa- 
tion (30+^ Myr) 2/ A more conservative, age range - 30 +}q° Myr 
- defined by the lower age limit from kinematics and the upper 
age limit from the empirical CMD. 

3.2. Comparison of k And b to reference objects 

We first used the photometry of k And b derived in Section 2 
to study the companion location in color-color (Figure 5) and 
color-magnitude diagrams (Figure 6). The new J, H, K s -band 
based colors of k Andromeda b are close to those of other late- 
M to mid-L class companions (US CO CTIO 108B, M9.5, Bejar 
et al. (2008); 2MASS J01033563-5515561(AB)b (Delorme et al. 

2013) ; GSC 06214-00210 b, L0, Ireland et al. (2011); AB Pic 
b, L0, Bonnefoy et al. 2013b; /3 Pictoris b, ~L2, Bonnefoy et al. 
(2013a); 1RXS J 160929. 1-2 10524b, L4, Lafreniere et al. (2008); 
CD-35 2722B, L4, Wahhaj et al. (2011)). These colors are also 
similar to those of young L0-L3 dwarfs (L r dwarfs; Kirkpatrick 
2005) listed by Faherty et al. (2013), and of field L dwarfs with 
availiable 2MASS photometry listed in Burgasser et al. (2006), 
Reid et al. (2008), and Dupuy & Liu (2012). The companion still 
lies close to GSC 06214-00210 b and to L4-L6 dwarfs in K s - L' 
vs J - K s and K s - L' vs H - K s diagrams. The companion col- 
ors are intermediate between those of 1RXS J 160929. 1-2 10524b 
and AB Pic b in these diagrams, k And falls intermediate be- 



H-K s 

Fig. 6. Position of k And b (magenta dot) in a color-magnitude 
diagram based on an improved near-infrared photometry. 


tween the sequence of field dwarfs and of the 1-3 Myr old com- 
panion GSC 06214-00210 b in a K s - L' vs L' - M' diagram. 
GSC 06214-00210 b is surrounded by circum(sub)stellar mate- 
rial, and shows signs of ongoing accretion (Bowler et al. 2011; 
Bailey et al. 2013). Therefore, the location of k And b in this 
diagram suggests that it has an intermediate age between GSC 
06214-00210 b and the population of field objects, or that it does 
not exhibit excess emission as strong as for GSC 06214-00210 
b. 

k And b falls on the sequence of L0-L4 field dwarfs in 
color-magnitude diagrams. The two L0 r objects of Faherty et al. 
(2013) with trigonometric parallaxes bracket the companion 
photometry. This is also self-consistent with the good match 
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Fig. 7. k And b spectral energy distribution (magenta dots) com- 
pared to the mean spectral energy distribution of young M8-L0 
dwarfs (Luhman et al. 2010) and of L0-L5 r dwarfs (Faherty et al. 
2013). 


to AB Pic b (classified as L0 following the good match with 
the young LOy dwarf 2MASS J01415823-4633574; Kirkpatrick 
et al. 2006; Bonnefoy et al. 2010). The companion falls above 
the 12^ Myr old exoplanet J3 Pictoris b, and at an intermedi- 
ate location between those of other known 30 Myr old compan- 
ions to 2MASS J01033563-55 15561 AB and HR 8799. This is 
consistent with k And being more massive than the exoplanets 
HR 8799bcde (<7-10-10-10 M Jup ; Marois et al. 2010; Currie 
et al. 2011; Sudol & Haghighipour 2012) and less massive than 
2MASS J01033563-5515561(AB)b (12-14 M Jup ; Delorme et al. 
2013), if we assume the system is a member of the Columba 
moving group (see Section 3.1). 

We finally analyze the companion SED from the available 
photometry, the corresponding filter pass bands, and a flux cal- 
ibrated spectrum of Vega (Bohlin 2007), as shown in Figure 7. 
We also built, and report on the figure, the SED of typical young 
M8-L0 dwarfs using the 2MASS and Spitzer colors 5 reported in 
Table 13 of Luhman et al. (2010) and normalized to the K s band 
flux of k And b (the error on the normalization factor were prop- 
agated). We reconstructed the SED of L r dwarfs taking the mean 
of the photometry reported in Table 3 of Faherty et al. (2013) for 
a given spectra type, and normalizing the flux to the one of k 
And b in the K s band . The SED of k And b falls between those 
of typical M9 and L0 dwarfs. It is incompatible with SEDs of L0 
and L3 r dwarfs. 

In summary, our empirical analysis suggests that k And 
b is a M9 to L3 dwarf. We use the bolometric correction 
BCk of 3.39+0.03 mag (Todorov et al. 2010) derived for the 
prototype L0 r dwarf 2MASS J01415823-4633574 (Kirkpatrick 
et al. 2006, and considering the BCk of the Taurus M9.5 
member KPNO-Tau 4 and of another young L0 field dwarf 
2MASS J0241 11510326587 for the error bar) in order to esti- 

5 We used the Gemini magnitude-flux converter for that purpose 
(http://www.gemini.edu/sciops/instruments/midir-resources/imaging- 
calibrations/fluxmagnitude-conversion). 


Table 3. Best- fit atmospheric parameters for k And b 


Atmospheric model 

Teff 

(K) 

logg 
(cm.s 2 ) 

R 

(Rjup) 

+ 

AMES-Dusty 

1900 

4.5 

1.25 

7.53 

AMES-Cond 

1700 

3.5 

1.47 

32.02 

BT-Settl 2010 

2000 

5.5 

1.18 

10.25 

BT-Cond 2012 

1800 

4.0 

1.34 

33.55 

BT-Dusty 2012 a 

1800 

4.5 

1.65 

14.10 

BT-Settl 2012 [M/H]=0.0 

1900 

4.0 

1.26 

12.05 

BT-Settl 2012 [M/H]=+0.5 

1900 

4.0 

1.24 

12.59 

DRIFT-P. [M/H]=0.0 

2000 

3.5 

1.13 

13.29 

DRIFT-P. [M/H]=+0.5 

1900 

3.5 

1.16 

12.35 

DRIFT-P. [M/H]=-0.5 

1900 

3.5 

1.16 

16.42 


Notes. (fl) Analysis limited to log g > 4.5. 


mate Logio(L/L©) = -3.76 ± 0.06 for k And b. This bolometric 
correction remains close to the one valid for M9-L8 field dwarfs 
(BC k = 3.19 - 3.33 mag; Liu et al. 2010). 

3.3. Atmospheric models 

We compared the SED of k And b to synthetic fluxes generated 
from seven atmospheric grids (AMES-DUSTY, AMES-COND, 
BT-DUSTY, BT-COND, BT-Settl 2010, BT-Settl 2012, DRIFT- 
PHOENIX) in order to determine the atmospheric parameters 
of the companion (T e ff, log g) and evaluate systematic errors 
on these parameters introduced by the models (Bonnefoy et al. 
2013b, submitted). The models and the fitting procedure are de- 
scribed in Bonnefoy et al. (2013a). We used the 2010 release 
of the BT-Settl models (BT-Settl 2010), despite the models do 
not incorporate the up-to-date physics. Indeed, Bonnefoy et al. 
(2013b) showed that BT-Settl 2010 synthetic spectra tend to bet- 
ter reproduce the near-infrared (1.1 -2.5 pm) spectra of young 
objects at the M-L transition, such as k And b. Results are re- 
ported in Table 3. The best fitted synthetic fluxes are displayed 
in Figure 8. x 2 maps of the fit (with 3 and 5 cr confidence levels 
overlaid) are shown in Figure 9. 

The models successfully reproduce the companion SED for 
T e ff = 1900^00 K. This corresponds roughly to the extreme val- 
ues determined for a 3cr confidence level using the most ad- 
vanced models BT-SETTL 12 (Figure 9). The temperature coin- 
cides with the latest range derived for the young L0 objects AB 
Pic b and 2MASS J01415823-4633574 from near-infrared spec- 
tra and SED fits (T eff = 1800^oo K; see Bonnefoy et al. 2013b). 
We also find that models with photospheric dust produce bet- 
ter fits to the data. These models provide a better constraint on 
T ef f, but not on the surface gravity (log g=4.5 ± 1.0). This was 
already the case for Pictoris b (Bonnefoy et al. 2013a). The 
high content in photospheric dust of k And b is also consistent 
with constraints on the spectral type derived in Section 3.2. To 
conclude, we note that the degenerate effects of metallicity, and 
surface gravity do not affect the temperature determination of 
the companion. 

The effective temperature and bolometric luminosity de- 
rived in Section 3.2 give a semi-empirical radius estimate of 
1 .2 ^qj Rj up for k And b. This radius is consistent with the value 
derived by adjusting synthetic fluxes to the observed companion 
SED exepted for the BT-DUSTY models, which temperature de- 
termination is likely biased by the limited coverage of the grid 
(Table 3). The radius matches predictions from “hot- start” evo- 
lutionary tracks corresponding to the companion luminosity and 
temperatures for ages between 30 and 250 Myr. 
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Fig. 8. Best fitted synthetic flux (horizontal bars) to the spectral-energy distribution of k And b (magenta dots) considering the 
DRIFT-PHOENIX models with three different metallicities (left panel), the BT-SETTL models (middle panel), the AMES-DUSTY, 
AMES-COND, and BT-COND models (right panel). The corresponding synthetic spectra are overlaid in each panel. 


Table 4. Mass of k And b prediced by “hot- start” evolutionary 
models 


Model 

age 

Mass from T eff 

Mass from L/L 0 


(Myr) 

(Mj up ) 

(Mjup) 

COND03 

COND03 

FM08 

SB12 

30-iq 

30 +1 ^° 

3^8 

30^8 
-10 

14 + “ 

14 + f 

> io 

> 10 

13 + 1 
13 + f 
> 10 
> 10 


3.4. The mass of k And b 

We compared T e ff and the luminosity estimates derived in pre- 
vious sections to predictions of evolutionary models in order to 
re-estimate the companion mass. We considered the two distinct 
age ranges determined in Section3.1. 

3.4.1. Classical hot- and cold-start models 

We first used the “hot-start” models of Baraffe et al. (2003) 
(hereafter COND03). We compare in Figure 10 masses predicted 
by these models for different companions, k And b mass esti- 
mates are in the brown-dwarf regime if the system is 150 Myr 
old. 

We used alternatively the “hot-start” models of Fortney et al. 
(2008) (FM08), and Spiegel & Burrows (2012) (SB 12). FM08 
and SB 12 models explore the impact of chemical enrichment re- 
lated to the formation process (lx and 5x solar for the FM08 
models, lx and 3x solar for the SB 12 models) on the object 
emergent flux and evolution. Results are reported in Table 4. 
These models predict masses above 10 M Jup (maximum mass 
covered by these models), k And b’s temperature and bolometric 
luminosity give the same estimates. 

We also considered the “cold-start” version of FM08 and 
SB 12 models. Model predictions do not extend to sufficient high 
masses (M > 10 M Jup ) to reproduce the luminosity and tempera- 
ture of k And b. 


3.4.2. Warm-start models 

The formation mechanism of k And b is not known (see section 
4) and, more importantly, the outcome, in terms of initial 
brightness, of the different formation scenarios cannot yet be 
predicted. We show in Figure 10 the impact of initial conditions 
with two cooling curves of Marleau & Cumming (2013). The 
“cold-start” curves correspond to cases with initial entropie 
of 9.5 Boltzmann units per baryon (k B /baryon). The “hot- start” 
cooling curves correspond to Sj n i t =13 k B /baryon for a mass of 
3 M Jup , and 14 k B /baryon for masses from 5 to 13.6 M Jup . It 
is then essential to take into consideration models with a wide 
range of possible entropies, including those of the “hot-start” 
and “cold-start” models. We examined the predictions of two 
sets of “warm- start” models (Spiegel & Burrows 2012; Marleau 
& Cumming 2013) for k And b for that reason. 

The “warm-start” models of SB 12 consider from 8 to 
13 k B /baryon, in 0.25 k B /baryon increments, and masses from 1 
to 15 Mj up . The models incorporate a deuterium-burning phase 
(Spiegel et al. 2011a, A. Burrows, priv. com.). Temperature and 
radii predictions were provided by the authors (D. Spiegel priv. 
com.) for this set of input parameters. We combined them to 
create bolometric luminosity predictions. The luminosity and 
effective temperature predicted by the models tend to increase 
with the object mass and initial entropy, and decrease with ages. 
Therefore, we determined the combination of initial entropies 
and masses corresponding to k And b’s measured luminosity and 
temperature assuming a system age of t = 30 Myr and propagat- 
ing the associated uncertainties. We show the results in Figure 
1 1 for an age of t - 30+^ Myr. The mass of k And b is greater 
or equal to 13 M Jup according to these models. Predicted masses 
from the luminosity agree with the ones derived from the tem- 
perature, althougth predictions do not extent to sufficient high 
masses to reproduce the upper limit on the estimated T e ff of k 
And b. The companion properties can also not be reproduced for 
ages of 150 Myr for the same reason. 


10 


M. Bonnefoy et al.: Characterization of the gaseous companion k Andromedae b 



l°9(g) (dex) 



109(9) (dex) 



109(9) (dex) 





109(9) (dex) 






Fig. 9. X 1 maps corresponding to the comparison of the spectral energy distribution of k And b to synthetic fluxes derived from 
atmospheric models for given log g and T e ff. Minima are indicated by magenta dots. We overlay contours corresponding to 3cr 
(green) and 5 cr (red) confidence levels. 
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Fig. 10. Evolution of the luminosity of gaseous objects pre- 
dicted by the COND models (black solid line), and Marleau 
& Cumming (2013) models with typical “hot-start” (light pink 
dashed curve; 3, 5, 10, 13.6 M Jup ), and “cold-start” initial condi- 
tions (dark blue dashed curve; 3, 5, 10, 13.6 M Jup ). We over- 
lay measured luminosity of young low mass companions. A 
more complete version of this figure can be found in Marleau 
& Cumming (2013). 


We also derived absolute flux predictions of SB 12 models 
for the given filter passbands and the four sets of boundary 
conditions (cloud-free models at solar metallicity - cfls, cloud- 
free models with three times the solar metallicity - cf3s, hybrid 
clouds at solar metallicity - hyls, hybrid clouds with three times 
the solar metallicity - hy3s) used for k And b following the same 
method as in Bonnefoy et al. (2013a). These synthetic fluxes 
were compared to the observed SED. The results are reported in 
Table 5. The comparison is biased by the limited mass coverage 
of the models. We note however, that solutions found within the 
models boundaries correspond to initial entropies intermediate 
between those of hot and cold-start models, placing the mass at 
the typical planet/brown-dwarf boundary (-13.6 Mj up Spiegel 
et al. 2011b; Molliere & Mordasini 2012; Bodenheimer et al. 
2013). These solutions correspond to T eff values that are in good 
agreement with those determined from the companion SED. 


Table 5. Best fit photometric predictions of the “warm-start” 
evolutionary models. Solutions found at the edges of the param- 
eter space (mass, Si n i t ) covered by the models are highlighted in 
italic. 


Atmospheric model 

Age 

(Myr) 

Mass 

(Mj up ) 

Sinit 

(k B /baryon) 

x 1 

Cloud free - lx solar 

20 

15 

9.75 

83.19 

Cloud free - 3x solar 

20 

15 

9.75 

57.77 

Hybrid cloud - lx solar 

20 

14 

9.75 

14.56 

Hybrid cloud - 3x solar 

20 

14 

9.75 

11.92 

Cloud free - lx solar 

30 

14 

9.75 

83.39 

Cloud free - 3x solar 

30 

14 

9.75 

58.32 

Hybrid cloud - lx solar 

30 

14 

10.25 

13.73 

Hybrid cloud - 3x solar 

30 

14 

10.00 

11.29 

Cloud free - lx solar 

50 

14 

10.25 

82.45 

Cloud free - 3x solar 

50 

14 

10.50 

57.66 

Hybrid cloud - lx solar 

50 

14 

13.00 

19.30 

Hybrid cloud - 3x solar 

50 

14 

13.00 

16.27 

Cloud free - lx solar 

150 

14 

13.00 

185.43 

Cloud free - 3x solar 

150 

14 

12.75 

175.71 

Hybrid cloud - lx solar 

150 

14 

13.00 

173.50 

Hybrid cloud - 3x solar 

150 

14 

13.00 

168.67 
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Fig. 11. Predictions of the “warm-start” evolutionary models of 
SB 12 for k And b for a system at t ~ 30^ Myr. The extreme val- 
ues of the companion age, T eff , and luminosities define a range 
of masses and initial entropies lying between the dashed and 
dotted-dashed curves. We also overlay the initial entropies con- 
sidered in “hot-start” (open circles) and “cold-start” (dots) mod- 
els of FM08 (based on Marley et al. 2007). 


i i i i i i i i i i i i i 


L/U 


Estimate 
Lower limit 
Upper limit 


i i i i i i i i i i i i i 


In comparison, the models of Marleau & Cumming (2013) 
have a much simpler outer boundary condition (hereafter 
MCI 3), using a grey, solar-metallicity atmosphere. We used 
them as they can be used to evaluate the impact of underly- 
ing hypotheses made in the models (e.g. atmosphere treatment, 
equation of state) on the derived joint mass and Si n i t values. We 
ran Markov Chain Monte Carlo simulations (MCMCs) in mass 
and initial entropy as in Marleau & Cumming (2013) with the 
related models to account for the uncertainties on the age, T eff , 


and luminosity of k And b. We assumed Gaussian distributions 
on L and T e ff . We took normal or lognormal errorbars for the two 
considered age ranges ( t = 30^ Myr and t = 30+ j^ 0 Myr), and 
chose flat priors in S i n i t and M. 

Fig. 12 displays the 68-, 95- and 99 % joint confidence re- 
gions from the MCMC runs for both age groups. Open and 
closed circles are as in Marley et al. (2007), show the approx- 
imate range of entropies spanned by hot and coldest starts, re- 
spectively, but shifted upwards by +0.38 k ^/ baryon to match the 
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Fig. 12. Joint posteriors on the mass and intial entropy of k And b 
obtained through MCMC runs with the models of Marleau & 
Cumming (2013), which include deuterium burning, from T q q 
(solid lines) and the bolometric luminosity (dashed lines) of k 
And b, assuming Gaussian errorbars. Entropy values need to be 
decreased by 0.45 & B /baryon when comparing with SB 12. The 
contours indicate the 68-, 95- and 99 % confidence levels (thick 
to thin lines). The colour (from dark to light) is proportional to 
the joint posterior probability at each M and S^t location in the 
r e ff- -based run. The age was taken as t - 30^ Myr (top panel) 
and t = 30+ J^ 0 Myr (bottom panel). 

luminosity in the models of MCI 3 (see therein). The results are 
consistent with those of Fig. 1 1 . 

Even for the low-age group, almost all solutions are in the 
mass regime where deuterium burning is important for the evo- 
lution of the object. A dramatic illustration of this lies in the 
solutions found at lower entropies. Whereas the models of SB 12 
allowed, for an age of 30+^ Myr, initial entropies down to only 
9.5 & B /b aryon to 1 cr in luminosity (see Fig. 11), the models 
of MCI 3 find that k And b could have formed with an entropy 
as low as « 8.8 & B /b aryon, correcting downward from Fig. 12 


IQ ” 2 



Time (yr) 

Fig. 13. Examples of possible cooling curves for k And b using 
log F/E© = -3.76 at 30 Myr, with the corresponding mass and 
initial entropy (to be decreased by 0.45 & B / baryon when com- 
paring with SB 12) of each curve indicated in the figure. For the 
lower entropies, we find that k And b could be presently un- 
dergoing a deuterium ‘flash’ (Salpeter 1992; Bodenheimer et al. 
2013, Marleau & Cumming, in prep.). 


for the entropy offset of 0.45 & B /baryon between the two mod- 
els 6 . The low-5 Mt solutions are possible only if the models in- 
clude a rise in the object’s luminosity due to deuterium burning 
(Molliere & Mordasini 2012; Bodenheimer et al. 2013, Marleau 
& Cumming in prep.). This is illustrated in Fig. 13, which shows 
cooling tracks for different (M, S^t) combinations which all 
reach logL/F© = -3.76 at 30 Myr. The lowest initial-entropy 
solutions undergo a ‘flash’, where the entropy in the object in- 
creases on a short timescale. The combinaison of the measured 
T e ff and luminosity (and therefore the object radius) can not help 
to discriminate these different possible cooling curves. Indeed, 
the high initial-entropy (“hot-start”) and the flashing cooling 
curves have a difference of some 25 K or 0.05 Rj in predicted 
r eff and radius only. 

The exact low-entropy solutions can depend on the details 
of deuterium burning, for instance on the initial deuterium con- 
tent or metallicity of the object (Spiegel et al. 2011b; Molliere 
& Mordasini 2012), but the main conclusion is a robust one: 
the combustion of deuterium may play a significant role in the 
cooling history of k And b, irrespective of the precise age of the 
system. We warn that while the lowest S Mt are comparable to the 
extrapolation of the coldest starts (Marley et al. 2007) to higher 
masses, what this implies about the formation mechanism is not 
clear given the major uncertainties about their outcome (see also 
the discussion of possible formation processes in Section 4). 


6 The implementation of the equation of state of Saumon et al. (1995) 
is slightly different in the two cases, with SB 12 using a simpler version 
of the Saumon et al. (1995) code that does not include a contribution 
from the proton spin in the partition function. See Marleau & Cumming 
(2013) for details. 
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Table 6. Properties of the k And system 


Parameter 

k Andromedae A 

k Andromedae b 

Ref 

d (pc) 

51.6 + 0.5 


2 

Age (Myr) 

so:™ 


1,3 

J (mag) 

4.26 ± 0.04 

15.86 + 0.21 

1 

H (mag) 

4.31+0.05 

14.95 + 0.13 

1 

K s (mag) 

4.32 + 0.05 

14.32 + 0.09 

1 

L’ (mag) 

4.32 + 0.05 

13.12 + 0.1 

1,3 

NBA. 05 (mag) 

4.32 + 0.05 

13.0 + 0.2 

1 

M’ (mag) 

4.30 ± 0.06 (a) 

13.3 + 0.3 

1 

mj(mag) 

0.70 + 0.06 

12.30 + 0.22 

1 

m H (mag) 

0.75 + 0.06 

11.39 + 0.15 

1 

m Ks (mag) 

0.76 + 0.06 

10.75 + 0.11 

1 

m L ' (mag) 

0.76 + 0.06 

9.56 + 0.11 

1 

m NB_4.05 (mag) 

0.76 + 0.06 

9.44 + 0.23 

1 

m M ' (mag) 

0.85 + 0.06 

9.75 + 0.31 

1 

Spectral type 

B9IVn 

M9-L3: 

1,4 

Teff (K) 

10900 + 300 

1900 +10 ° 
1 VUU_200 

1 


11400+ 100 

5 


10700 + 300 


6 

log g (dex) 

3.78 + 0.08 (fc) 

4.5+ 1.0 

1 


4.10 + 0.03 


5 


3.87 + 0.13 


6 

M/H (dex) 

-0.36 + 0.09 


1,5 


-0.32 + 0.15 


6 

logio(L/L Q ) 

1.83 + 0.04 (c) 

-3.76 + 0.06 

1 

Mass (M©) 
Mass (M©) 

2.6 + 0.2 (J) 

0 01 'U 0 - 022 
u.uij_oooi 

> 0.011 (7) 

1 

1 


Notes. [1] - this work, [2] - van Leeuwen (2007), [3] Carson et al. 
(2013), [4] Garrison & Gray (1994), [5] Fitzpatrick & Massa (2005), 
[6] Wu et al. (2011). 

Estimated from the mean K’-M’ colors of B9 stars. 

The polar surface gravity. Estimated by correcting the measured 
log(g) for the rapid rotation of the star using the method of Huang & 
Gies (2006). 

(c) Found ajusting BT-Settl spectra with T e ff =10 600, 10 800, and 11 
200 K, and log g=3.5 on the 0.365-22.1 /um flux distribution of the 
source. 

{d) Estimated comparing the measured T eff and bolometric luminosity 
of k And A to evolutionary tracks of Ekstrom et al. (2012), with and 
without rotation, for ages of 20 to 150 Myr. 

(e) Using COND “hot-start” evolutionary models for the most concer- 
vative age range (Baraffe et al. 2003). 

^ Using the “warm- start” models. 


4. Discussion 

The characterization of k And further illustrate the challenge of 
determining accurate masses for companions due to uncertain- 
ties on the age-dating methods and evolutionary tracks. 

The de-projected (and projected) separation of k And b 
(61^0 AU, Carson et al. 2013) is compatible with the size of 
primordial (e.g. Mannings & Sargent 1997; Panic & Hogerheijde 
2009; Guilloteau et al. 2013, and ref therein) and debris disks 
(Booth et al. 2013, and ref therein) surrounding stars in the same 
mass range as k And A, some of which show structures sug- 
gesting a clear signpost for planets (e.g. HR4796, Schneider 
et al. 1999; Lagrange et al. 2012b). The companion location also 
fits well with the extent and the location of cavities/gaps/spirals 
in young structured (or transition) disks (Andrews et al. 2011; 
Grady et al. 2013; Quanz et al. 2013b) discovered around Herbig 
Ae stars, k And b’s separation is close to that of the candidate 
substellar embryo (Quanz et al. 2013a) around the 2.4 M© star 
HD 100546, which mass might extend inside the “brown-dwarf’ 


regime. It is in addition intermediate between the two outermost 
planets orbiting HR 8799, and nearly identical to the proba- 
ble exoplanet HD 95086b (Rameau et al. 2013a), all orbiting 
intermediate-mass stars. Therefore, despite the large uncertain- 
ties on the mass of k And b and the lack of disk excess emission 
around k And A (Section 3.1), we should still consider that the 
companion could have formed within a disk. 

The growing population of massive gaseous companions dis- 
overed on short-period orbits around massive stars (e.g. Galland 
et al. 2006; Deleuil et al. 2008; Hartmann et al. 2010; Lee 
et al. 2011; Bouchy et al. 2011, and ref therein), the discov- 
ery of unusually dense substellar companions (e.g. CoRoT-20b 
and HAT-P-20b Bakos et al. 2011; Deleuil et al. 2012), and 
recent simulations (Mordasini et al. 2009b; Bonnefoy et al. 
2013a), suggest that the core- accretion mechanism might still 
work well inside the brown-dwarf mass regime 7 . Kennedy & 
Kenyon (2008) shows that 10 M® cores can form in less than 
1 Myr from ~3 to ~23 AU around 2.4-2. 8 M® stars. Rafikov 
(2011) also propose that core- accretion could operate out to 40- 
50 AU. A formation closer to the snow line associated to dy- 
namical scattering or outward migration could explain k And b 
(Paardekooper & Papaloizou 2008; Peplinski et al. 2008; Crida 
et al. 2009b; Paardekooper et al. 2010; Chatterjee et al. 2011; 
Lin & Papaloizou 2012). Dynamical scattering would neverthe- 
less require additional, but yet undetected massive companions. 
Therefore, it is less suitable if the system is significantly older 
than 30 Myr and/or if the initial entropy of k And b is close 
to “cold-start” conditions. The proposed modifications on the 
way solids are accreted to form a core (Ormel & Klahr 2010; 
Lambrechts & Johansen 2012; Morbidelli & Nesvorny 2012) 
might still facilitate the formation of k And b by nucleated in- 
stability closer to its present separation. 

We also applied a dedicated formation model to explore the 
possibility that k And b could originate from a disk-instability 
(Klahr et al. in prep; see also Janson et al. 2011; Rameau et al. 
2013b, for a description of the models). These models can pre- 
dict the range of semi-major axes and companion masses orig- 
inating from clumps that 1/ formed in Toomre-unstable disks 
(Toomre 1981) and, 2/ cooled down more rapidly than the local 
Keplerian timescale. The models require as input the initial lumi- 
nosity of the star and the system metallicity. The initial luminos- 
ity was estimated by inputting the temperature and luminosity of 
the star (see Table 6) into evolutionary models of Ekstrom et al. 
(2012), while considering system ages of 20 and 150 Myr. We 
generated disk models with solar and under- solar abundances in 
order to reflect the possible metallicities of the star (Note that the 
star probably has a solar metallicity, see Carson et al. 2013). The 
two metallicities do not change significantly the predictions (see 
Figure 14). In contrast with Pictoris b (Bonnefoy et al. 2013a; 
Currie et al. 2013), or HR 8799bcde (Janson et al. 2011), k And b 
has an estimated mass and semi-major axis (and projected sep- 
arations) compatible with a formation by disk-instability close 
to its present location. The masses of k And b fall at the edge 
of the allowed range set by the Toomre criterium if the system 
age is 20 Myr. Adopting the lower limit on the estimate of the 
initial luminosity (L/L® =1.55) shifts the contrain associated 
to the Toomre criterium to lower mass values, therefore solv- 
ing the issue. This also highlight the fact that these models pre- 
dictions remain highly dependent of the correct determination 


7 We note that disk-instability associated to “tidal-downsizing” might 
also explain the observations (Boley et al. 2010, 2011; Nayakshin 
2010b, a, 2011). But the models hasn’t been extensively tested against 
observations yet (Forgan & Rice 2013). 
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Fig. 14. Predicted masses of fragments caused by disk instabili- 
ties as a function of semi-major axis for the massive star k And. 
Formation of planets by disk instability is possible when the 
Toomre criterium is exceeded (solid lines) and when the disk 
can cool sufficiently rapidly (long dashed lines). The allowed 
formation space (in yellow) lies in-between these two criteria. 
Companions can not be formed in situ at separations shorter than 
the ones indicated by the dotted lines. We superimpose the semi- 
major axis and mass estimates of k And b (magenta dot) esti- 
mated from “hot-start” models for Columba ages (30+^ Myr; 
dark grey shaded area) and considering an upper-limit corre- 
sponding to the highest age for the Pleiades cluster (3 CC|q° Myr; 
light grey shaded area). We also report the lower limit on the 
companion mass inferred from “warm-start” models (grey ar- 
rows, see Section 3.4.2). We show the corresponding primordial 
disk-mass reservoir (10, 20, 50%) at a given separation needed 
to form given clumps/companions. The Toomre parameter is in- 
dependent of stellar metallicity. 


of these initial luminosities. Recent alternative disk-instability 
models support the idea that fragments tend to migrate inward on 
extremely short timescales (e.g. Michael et al. 2011; Baruteau 
et al. 2011; Zhu et al. 2012; Vorobyov 2013, and ref therein), 
thus arguing against an in- situ formation by disk-instability for 
k And b. Nevertheless, our model also show the companion - if 
massive enougth - might still be formed up to 220 AU. Also, mi- 
gration processes are still far from being understood (Mordasini 
et al. 2009b). 

We applied the same models to the current population 
of directly-imaged brown-dwarfs companions orbiting young 
2.2-2.5 M q stars (HR 7329 B, HD 1160B, HIP 78530B; see 
Appendix C) in an attempt to understand how k And b relates to 
these objects. HR 7329 B has estimated separation and masses 
(using “hot-start” models) compatible with a formation by disk- 
instability. This companion and k And b fall in a similar area 
in the diagrams (Figures 14 and C.l) where the primordial disk 
needed to form clumps with the estimated masses of the compan- 
ions is relatively light (~10% M*). Such disk masses are com- 
patible with the observations of primordial disks (e.g. Beckwith 


et al. 1990; Gaczkowski et al. 2013). HR 7329 B is the only 
companion of the sample known to orbit a star with a debris 
disk (see Su et al. 2006; Chen et al. 2012, for HD 1160 and 
HIP 78530). Smith et al. (2009) derive an outer radius of 24 
AU, which might result from truncation/disruption (Rodriguez 
& Zuckerman 2012). HD 1160 B is part of a multiple system 
with a 0.22 Ms un companion HD 1160 C imaged at a projected 
separation of 533 AU (2.2”) from the star, therefore possibly 
placed on a concentric orbit. Models indicate that HD 1160 B 
would need a more massive (~0.2 M star ) disk to form than k And 
b. Nevertheless, the comparison of HD 1160 B properties to the 
disk-instability models is more subject to caution since the for- 
mation of companions in such high-multiplicity systems is prob- 
ably far more complex. To conclude, the model predict that the 
extremely wide companion HIP 78530 B Lafreniere et al. (201 1) 
can not have formed in- situ, contrary to the three other studied 
systems. Therefore, these comparisons would now deserve being 
repeated with other models, and on a larger sample of objects. 

5. Conclusions 

We present the first deep-imaging observations of the k And sys- 
tem at 4.05 (NB_4.05) and 4.78 pm (M'). We retrieved the com- 
panion at these wavelengths and estimate NB_4.05 = 13.0 ± 0.2 
and M' = 13.3 + 0.3 mag. We also obtained new photometry of k 
And A in the J, H, and K s bands, which we use to re-evaluate the 
1 .1-2.5 pm photometry of the companion. The resultant photom- 
etry indicates that the companion is a late-M or early-L dwarf 
with a bolometric luminosity of Logio(L/L©) = -3.76 + 0.06. 
We obtained a low-resolution set of observations accross the 1- 
5 pm spectral energy distribution of k And b and compared it 
to predictions from atmospheric models exploring the limiting 
and intermediate cases of dust-formation. All models converge 
toward a T e ff = 1900^oo K for k And b. Models with dust in 
the photosphere of the object better reproduces k And b spectral 
energy distribution. The models do not enable to constrain the 
companion’s surface gravity. 

The luminosity and temperature were then used as input 
of evolutionary models accounting for a wide range of initial 
conditions. We re-estimated and considered for that purpose a 
more conservative age range than in the discovery paper (30+jg° 
Myr). “Hot- start” models constrain the mass of k And b to 
12-39 Mj up . Conversely, “Warm-start” models computed for 
initial entropies from 8 to 13 k B /baryon provide a lower limit 
of 11 M Jup . Therefore a substantial fraction of the allowed mass 
range of k And b is in the typical brown dwarf regime (> 13.6 
Mj up ). The latest “Warm-start” models reveal in addition that k 
Andromedae b could be undergoing a deuterium flash. This flash 
is expected to play a significant role in the cooling history of the 
companion. It poses a serious challenge to the companion mass 
determination, irrespective of the uncertainties associated to the 
system age. 

The formation models we used indicate that, for a large frac- 
tion of plausible masses, and given the separation of the com- 
panion, k And b might have formed close to its present loca- 
tion by disk-instability. We apply the same models to the cur- 
rent sample of low-mass companions on wide orbits (» 10 AU) 
around stars in the same masse range. These models suggests 
that some, but not all, of these objects could have also formed 
by disk-instability. 

Konopacky et al. (2013) has recently attempted a determina- 
tion of the C/O ratio in the atmosphere of the young wide-orbit 
gas giant planet HR 8799c. Applying the same method to k And 
b could help to clarify the companion’s formation scenario in 
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the near future. Additional measurements (radial velocity, non- 
redundant masking) would also be of value to set constraints 
on possible dynamical perturbers at smaller orbital separations, 
which could account for the current location of k And b. 
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Appendix A: k And A color-magnitude diagram 
samples 

We describe here each of the empirical sub-samples that were 
used to construct Figure 3: 

- Scorpious-Centaurus - For this sample, we used the 
list of Chen et al. (2012). We removed known spectro- 
scopic, eclipsing, and sub-arcsecond resolved binaries and 
Herbig AeBe stars to minimize photometric scatter in the 
CMD. Photometry and distances were collected from the 
Hipparcos catalog (Perryman et al. 1997; van Leeuwen 
2007). Photometry was corrected using the individual extinc- 
tions listed in Chen et al. (2012). 

- IC 2391 - This sample is comprised of early-type stars se- 
lected from the membership list of Perry & Bond (1969). 
Distances and photometry are from Hipparcos. The photom- 
etry was corrected using an average cluster reddening of 0.01 
(Patten & Simon 1996) 

- The Pleiades- We drew early-type members from the list of 
Stauffer et al. (2007) and used individual distances and pho- 
tometry from Hipparcos. Individual reddening values from 
Breger (1986) were used to correct the photometry. 

- Ursa Majoris moving group - We chose A-type stars pro- 
posed to be members of the group nucleus by King et al. 
(2003). Photometry calculated in this reference was also 
adopted. We applied no correction for reddening. 

- Young Moving Groups - We compiled proposed A and 
B-type members of the AB Doradus, Tucana/Horologium, 
Columba, and j3 Pictoris young kinematics groups from Malo 
et al. (2013, and references therein). For these stars, when 
available, we compiled the mean photometry from the cata- 
logs of Mermilliod (Mermilliod & Mermilliod 1994, includ- 
ing k And A). When photometry was not available in these 
catalogs it was taken from Hipparcos. In the construction of 
the CMD, we used Hipparcos distances (van Leeuwen 2007) 
and applied no reddening correction. 

Appendix B: New determination of the atmospheric 
parameters of k And A 

Our APO-ARCES spectrum of k And A resembles those of other 
of luminosity class IV and V B-type stars. Our spectral anal- 
ysis followed the procedure developed by Marsh Boyer et al. 
(2012) which uses model spectral template fitting to estimate 
the atmospheric parameters of the star. Their method uses four 
lines, Hy 44340 A, He I 44387 A, 44471 A, 44713 A, and Mg 
II 44481 A to determine vsin /, T e ff, and log g (see Figure B.l). 
Since the lines that are sensitive to vsin i (He I and Mg II) were 
weak in our spectrum, we maintained a fixed vsin / (vsin i =190 
km s" 1 , Glebocki et al. 2000; Fitzpatrick & Massa 2005) with a 
conservative uncertainty (20 km s -1 ) during the analysis of the 
Hy line, which is sensitive to temperature and gravity. This anal- 
ysis found two solutions: one where the model spectrum was a 
better fit to the line wings and core, but was slightly off in the line 
width, and another where the best fit model better reproduces the 
line width but gave a poorer fit to the wings. We combined the re- 
sults from these two solutions to estimate T e ff = 10900 ± 300 K, 
log g = 3.50 ± 0.08 dex where the larger uncertainties reflect the 
combined results. As described in Section 3, we applied a cor- 
rection to the measured surface gravity for the rapid rotation of 
the star using the method of Huang & Gies (2006). This method 
uses models to estimate the surface gravity at the pole of the 
star which should remain relatively unaffected by rapid rotation 



A (A) 



A (A) 


Fig. B.l. Hydrogen H r line of k And A (solid lines) fitted by 
atmospheric models (dotted lines). Residuals from the fit are 
shown in the upper part of the figures. Top: Solution where the 
line wings and the core are well fitted, but where the fit degrades 
in the line breadth. Bottom: Solution for which the line breadth 
is better fitted. These solutions were combined to derive the new 
estimate of the temperature and the surface gravity of the star. 

an give a better indication of the stars true evolutionary state. 
We find log g po i = 3.78 ± 0.08 dex for k And A. Our estimated 
T e ff is generally in agreement with those determined using other 
methods in the literature (e.g. Wu et al. 2011). However, the sur- 
face gravity is lower than previous estimates, which have not 
been corrected for the star’s rotation (see Table 6). The spread 
in estimated surface gravities illustrates the inherent challenges 
in accurate atmospheric parameter determination for early-type 
stars. 

Appendix C: Disk-instability models for HR 7329B, 
HD1160B, and HIP 78530B. 

We applied the disk-instability models of H. Klahr to the rare 
brown dwarfs companions identified around young (age <100 
Myr) early- A/late-B type stars. The models require as input the 
metallicity and initial luminosity of the star (the latter roughly 
scales with the stellar mass). We retrieve the zero-age main se- 
quence luminosity of the stars by inputting the present effective 
temperature or mass estimates of the stars and the known age of 
the system as inputs of evolutionary tracks Ekstrom et al. (2012). 
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The properties of the systems are summarized in Table C.l. We 
considered a solar metallicity for HR 7329B as the present mea- 
surements are roughly solar (a variation of +0.17 dex does not 
affect the cooling time-scale significantly; Saffe et al. 2008). We 
also find a small sensitivity of the modelsimulation of predic- 
tions to the choice of the metallicity for HD 1160 (Fe/H=0.0 and 
-0.3 considered here; see section 3.4 of Nielsen et al. 2012). 
Finally, we assumed HIP 78530B has a roughly solar metallic- 
ity. This is likely to be the case given recent measurements on 
lower mass stars of the associations (Viana Almeida et al. 2009). 

We deprojected the observed separation of HIP 78530B, HD 
1160B and C following the values from the Monte Carlo sim- 
ulation of Allers et al. (2009). We use the values reported in 
Neuhauser et al. (2011) for HR 7329 B. The properties of the 
three companions are compared to models predictions in Figure 
C.l. 
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Table C.l. Young 2. 2-2. 5 M© stars with brown-dwarf (or low-mass) companions on wide orbits. 
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age 

(Myr) 

FeH - star 
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Spectral 
type A 

Spectral 
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M8 + 1 
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23 + 3 
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6,7 


Notes. [1] - Neuhauser et al. (2011), [2] - Tetzlaff et al. (2011), [3] - Saffe et al. (2008), [4] - Nielsen et al. (2012), [5] - van Leeuwen (2007), [6] - 
This work, [7] - Lafreniere et al. (201 1). 

(fl) Estimated from “hot-start” models. 

(b) De-projected semi-major axis. 
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Fig. C.l. Same as Figure 14, but for the brown-dwarf companions HR 7329B (left, red cross), HD 1160B (middle, green triangle), 
and HIP 78530B (right, pink hourglass). We added the NaCo detection limit obtained by Rameau et al. (2013b) for HR7329 B. 
Model predictions do not extend beyond 1000 AU. 
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